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ABSTRACT 

Galaxy Zoo is the first study of nearby galaxies that contains reliable infor- 
mation about the spiral sense of rotation of galaxy arms for a sizeable number 
of galaxies. We measure the correlation function of spin chirality (the sense in 
which galaxies appear to be spinning) of face-on spiral galaxies in angular, real 
and projected spaces. Our results indicate a hint of positive correlation at sepa- 
rations less than ^0.5 Mpc at a statistical significance of 2-3 a. This is the first 
experimental evidence for chiral correlation of spins. Within tidal torque theory 
it indicates that the inertia tensors of nearby galaxies are correlated. This is 
complementary to the studies of nearby spin axis correlations that probe the 
correlations of the tidal field. Theoretical interpretation is made difficult by the 
small distances at which the correlations are detected, implying that substruc- 
ture might play a significant role, and our necessary selection of face-on spiral 
galaxies, rather than a general volume-limited sample. 
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1 INTRODUCTION 

Understanding the creation and evolution of the an- 
gular momentum of dark matter halos and galaxies 
is a crucial building block of a comprehensive theory 
of galaxy evolution. iHovld (|l949t ) was first to propose 
that the galaxy spin can be ascribed to the gravi- 
tational coupling with the surrounding galaxies. This 
idea h as b een forma lised and extended i n subsequent 
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the modern theory of the evolution of galaxy spin, kn own 
as the tidal torque theory (TTT; see ISchaefeiJ (|2008l ) for 
a review). This theory asserts that protohalos acquire 
most of their angular momentum in the early stages 
of their formation, from the lowest non- vanishing con- 
tribution from the linear Lagrangian theory, that is, a 
coupling of the quadrupole of the local mass distribu- 
tion to the external gravitational shear. Compared to 
A^-body simulations, theory produces qualitatively cor- 
rect results, although there are still significant discrepan- 
cies at a more quantitative level. Moreover, it seems that 
at present there are no clear theoretical directions for 
improving analytical models (jBarnes fc Efstathiou, 19871 : 
iPorciani eraLll2002al : iBailin fcSteinmet jbOOSl ). 

On the observational side, most of the work has been 
done using spiral galaxies. These are characterised by a 
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rotating disk of baryonic matter. The line perpendicular 
to the plane of the disk determines the axis of rotation, 
while the spiral arms in most galaxies encode the sense of 
rotation, i.e the difference between left-hand screw and 
right-hand screw sense of rotation. For spiral galaxies 
seen in projection, one can measure the observed galaxy 
ellipticities , which constrain the axis of the galaxy spin 
fPen ct al. 2000'; Lee fc "Penll2002l : iLee fc Erdogdull2007l : 
[Trujillo ct al. 2006i). This axis is known to within two-fold 
degeneracy associated with the tilt of the galactic plane 
with respect to the plane of the sky. Since the vector can 
point in two directions on the same axis, the elliptici- 
ties constrain the spin vector within a total of four-fold 
degeneracy. Note that chiral information, viz. informa- 
tion about the actual directions of the spin vectors as 
opposed to spin axis, is completely absent in the study of 
galactic ellipticities. However, this information contains 
important clues about the details of the emergence of the 
spin. As we will explain later in the text, the detection of 
chiral correlation function implies that the local inertia 
tensors must be correlated. This lends experimental sup- 
port to the theoretical expectations that the inertia and 
gravit ational shear tensor are correlated (jPorciani et al.l 
l2002bl ). 

We now have a unique tool to study the chiral prop- 
erties of galaxy spins. Through an online project called 
Galaxy Zoc0 ( Lintott et al. 2008), members of the public 
have visually classified the morphologies and spin orienta- 
tions for the entire spectroscopic sampl e of the Sloa.n Dig - 



ital Sky Survey (SDSS from now on) ([York et al" 



Data Release 6 (DR6) IjAdelman-McCarthv et al 



■20001 ) 
20081 ). 



The data an d its r eduction is extensively discussed in 
iLintott et al] l|2008l) . 

Spiral galaxies in the Galaxy Zoo sample are clas- 
sified as clockwise, anti-clockwise or edge on. The spin 
direction convention used here is such that clockwise 
and anti-clockwise rotations correspond to galaxies whose 
arms are r otating in the sen se of the letters Z and S re- 
spectively (jSugai fc IvdllQQSh . For each face-on galaxy we 
thus receive one bit of information corresponding to the 
sign of the galaxy spin vector projected along the line of 
sight. It is important to note that this information is inde- 
pendent of the tilt of the plane of the galaxy. We will refer 
to this one-bit information simply as galaxy spin. By the 
galaxy spin vector we mean the unit vector that defines 
the apparent spin of the galaxy: it is perpendicular to the 
disk plane and points in the direction the right turn screw 
would move if turned following the spiral arms inwards. 
This quantity is strongly correlated with the real angular 
momentum of the gas. The correlation, however, is not 
perfect and observations show that the angular momen- 
tum vector of the gas points in the opposite direction i n 
about four percent of systems (jPasha fc Smirnovlll982l ). 
In turn, there are theoretical expectations that there is 
a strong, but not perfect correlation between the angu- 
lar momentum vector of g as and that of the dark ma tter 
halos hosting the galaxy (jvan den Bosch eraLll2002l ). A 
detection of correlation in the galaxy spins would there- 
fore imply a correlation in the dark matter spin vectors. 
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Conversely, a non-detection of the spin correlation can 
be used to put upper limits on the correlation between 
angular momentum vectors of dark matter halos. 

This paper is structured as follows. In Section [2] we 
shortly review the tidal torque theory and its main re- 
sults. Section [3] will connect the correlation function rj to 
an observable correlation function of spins c, while the 
Section [4] will introduce our data and measurement tech- 
nique. We present our results and discuss systematics in 
Section [5l Finally, we discuss our results and conclude in 
the last Section |6l 



2 TIDAL TORQUE THEORY 

The tidal torque theory derives the following expres- 
sion for the angular momentum from the 1 st order lin- 
ear poturbation theory in Lagrangian space l|Whitell 19841 : 
iCatelan fc Th cuns 1996): 



Li(t) = a^{t)D€ijkTjiIik, 



(1) 



where a is the scale factor of the Universe, D is the growth 
factor and tijk is the Levi-Civita symbol. The local inter- 
tia tensor lij of the protohalo (the mass that will later 
form the dark matter halo) in Lagrangian space is given 
by 



3o / Qi 

Jv 



9jd q, 



(2) 



where qi are the Lagrangian coordinates around the cen- 
tre of mass of the halo and po is the mean density. The 
local shear tensor Tij is defined by 



didj(l>{ci) 



(3) 



where (p is the gravitational potential. In other words, 
the TTT requires two components: a non-vanishing 
quadrupole distribution of mass in the halo to be spun 
up and the cosmological tidal field. In principle, it sounds 
plausible to assume that while tidal fields between neigh- 
boring protohalos are correlated, since they are coming 
from the large scale modes, the local quadrupole moments 
of mass distribution are sourced due to random distribu- 
tion of the local inhomogeneities and should therefore be 
random. This assumption of a statistical isotropy of the 
inertia tensor, gives the following a natz for the angular 
moment correlator (|Pen et al.ll2000l ): 



—5ij + c — — TikTk- 



(4) 



where c controls the level of randomization of axial pref- 
erence due to non-linear and stochastic effects. In fact, 
this ansatz has been shown to satisfactorily explain the 
inclinations o f axes of spiral gala xies in vicinity of voids 
with c ~ 0.7 (|Truiillo et al.ll2006t ). 

Some further algebra gives the probability distribu- 
tion funct ion for spins s = L, usually assumed to be 
Gaussian l|Pen et al.ll2000l ): 



P(s|T) 



IQI 



-1/2 



4-K 



exp ■ Q ^ ■ ^) 



(5) 



Using this expression it is therefore possible to calculate 
various correlators of s if correlators of T are known. 
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We will now consider correlation functions. The 
most general form of possible correlation statistics of 
the galaxy spins consistent with the homogeneity and 
is otropy is the spin correlation tensor defined by (see e. 
g. iGroth et aLlll989l ) 

Eij{r) = (si(x)sj(x + r)> = 

[n(r)-E(r)]nf,+E(r)5,,. (6) 

Functions 11 and E are parallel and perpendicular corre- 
lation function. Following Porciani ot al. (2002a) we will 
be dealing exclusively with the "dot product" correlation 
functions given by 

J7(r)=H,, = (s(x)-s(x + r)), (7) 

and 

.,2(r) = ((s(x).s(x + r))^)-l/3. (8) 

It is important to realize what these two quantities mea- 
sure. The first one measures if the angular momentum 
vectors are correlated, while the second one measures if 
the axes of angular momentum vectors are correlated. 
Note that is perfectly possible to have ri2{r) > while 
ri{r) = 0. For example, if all spins were aligned along the 
z axis, but with an orientation that is chosen at random 
from z and —z, r]2 = 1, while 77 = 0. In fact, following the 
ansatz of Equations (|4} and ([5]) results in vanishing r](r) 
and a fin ite ri2(r) (analytic al expression for which can be 
found in iLee fc PenI (|200ll )). This is trivially seen from 
Equation ([5]), since P(s|T) — P(— s|T) and is a direct 
consequence of the assumption of isotropy of local mo- 
ments of inertia. To put it simply, for a fixed tidal field, 
averaging over possible realisations of the inertia tensor 
will, in general produce a preferred axis (determined by 
eigenvectors of T), but not a preferred direction. This 
is due to the fact that for every inertia tensor I that 
produces a final angular momentum L, an equally likely 
mirror image inertia tensor —I will produce an equal and 
opposite angular moment — L. Therefore, the common 
assumption that the local moments of inertia are ran- 
dom and uncorrelated will, in general, produce a non- 
vanishing axis correlation, but a vanishing correlation of 
the actual spin vectors. The important corollary is, that a 
detection of chiral correlations in the galaxy spins would 
directly indicate that the moments of inertia are non- 
random. 

3 CORRELATIONS OF PROJECTED SPINS 

3.1 Small scale correlation function of spins 

Our life as observers is complicated by the fact that the 
sense in which spiral arms wind in a projected image of 
a spiral galaxy measures the sign of the spin vector pro- 
jected along the line of sight rather than the spin vector 
itself. We will therefore consider the correlation function 
of s = sgn(s • z): 

c(r) = (S(x)S(x + r))) , (9) 

where we assumed that the radial vectors to positions x 
and x -I- r are parallel, i.e. the flat-sky approximation. 



Note that this only requires pairs of galaxies to be close 
enough so that the flat-sky approximation holds, rather 
than an entire survey occupying a small portion of the 
sky. For the time being we also neglect the difference 
between the dark-matter angular momentum and the gas 
angular momentum. 

To proceed, we note that r](r) is determined entirely 
by the one-point distribution function P[^\r) for cosine 
angle ^ = cos between the two spin vectors: 

nir) = {s,(x)s,(x + r)) = (m> = j f^P{f^\r) dfi (10) 

For a given /i, one will observe two galaxies with 
the same orientations of spins with the probability (see 
beginning of the Appendix 

p+i{e) = i-e/7v (11) 

and with different orientations of spins with the proba- 
bility 

P-i(6l) = l-P+i(6() = e/7r (12) 

The correlation function of the spin signs is then 
given by 

cir) = dMP(Mir) (P+i(e) - P-i{e)) , (13) 

where normalisation TV is in this case trivially given by: 

J dfiP{fi\r) (P+i (e) + P_i (61)) = 1 (14) 

When there are no correlations, P(/i) d/i = l/2d/i 
and both Tj{r) and c(r) = 0. When correlations exist, we 
must specify a one point probability distribution function 
for P{8\r). We assume the following form: 

P(6l|r) = sin(6') (1 + e(r)cos(6')) (15) 

Using this form, one obtains 

V{r) = ie(r) (16) 

c(r) = ie(r) (17) 

and so 

Vir) = ^c(r) (18) 

Equation (|18p hinges on the particular form for P{6\r) 
that we chose. In practice, different forms generically give 
the results that rj{r) — qc{r) with q typically between 1 
and 3/2. 

In reality, however, we measure the correlation func- 
tion only for galaxies that are sufficiently away from 
the edge-on orientation not to be classified as a face-on 
galaxy. For simplicity, let us assume that our sample con- 
tains only galaxies, whose spin vector satisfies 

s(x) ■ 2 > cos a. (19) 

In other words, galaxies that are inclined with an angle 
greater than a with respect to the line of sight are as- 
sumed to have been classified as edge-on. What is the 
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functional form for P+i{9) in this case? In Appendix lAl 
we show that for a > 7r/4 



7(61) e<iT-2a- 
f{9) ■K-2a<e <2a; 







and 



2a < 6 



< TT — 2q: 



(20) 







fi-Tv-e) n-2a<e < 2a- 



/(tt -e) 2a < I 



(21) 



where 



f{9) = 1 - 2cos(a) 



1 -1 / cos^ 
cos 

TT \ I 

2 COS a _i/cosq(cos6 
H cos ' 



sm U sm a 



(22) 



The "lost" probability, i.e. 1 — P+i — P_i corresponds to 
geometries that are not detected and in general result in 
A'^ < 1. Numerically integrating Equation (|13p . we can 
obtain a relation between rj and the measured c™°^''(r). 

3.2 Connection to gas angular momentum 

As discussed in IPasha fc Smirnovl (| 19821 ). a fraction / = 
0.04 of galaxies has gas angular momentum that is point- 
ing in the opposite direction to the apparent galaxy spin 
inferred from orientation of spiral arms. If we momen- 
tarily distinguish between the actual gas spin correlation 
function and apparent gas spin correlation function, we 
can write 

77appa.ont(r) = ((l-/)^ + /)r,(r)-2/(l-/)r?(r), (23) 

since correlation function of spins will receive a negative 
contribution if exactly one spin (but not both) was ran- 
domly reversed. This simplifies to 

77apparont (r) = ^ ^/ - r?(r). (24) 

This has the expected properties. The spin correlation 
function will become zero if exactly half the spin vectors 
are reversed, effectively randomizing them and exactly 
following the primary correlation function if all or none 
spins are reversed. For / = 0.96 one gets that ri{r) ~ 

1 .2?^ apparent (^) • 



4 DATA & METHOD 
4.1 Data 

T he basic data reduction i s described in gr eat detail 
in iLintott et al.l (|2008) and iLand etall (|200d). We wiU 
briefly summarise the data reduction in the following 
paragraph, but the reader is invited to read the above 
papers if interested in the details of the primary data 
reduction. 



In the Galaxy Zoo project, a sample 893,212 galax- 
ies were visually classified by about 90,000 users. The 
sample was selected to be sources that were targeted for 
SDSS spectroscopy, that is extended sources with Pet- 
rosian magnitude r < 17.77. Additionally, we included 
objects that were not originally targeted as such, but were 
observed to be galaxies once their spectrum was taken. 
Where spectroscopic redshifts are available, we find that 
they have the mean redshift of z = 0.14 and the ob- 
jects with the highest redshift reach z ~ 0.5. The galax- 
ies thus probe our local universe at cosmological scales. 
Each object has been classified about 40 times from a 
simplified scheme of 6 possible classifications; an ellip- 
tical, a clockwise spiral galaxy, an anti-clockwise spiral 
galaxy, an edge-on spiral galaxy, a star / unknown ob- 
ject, a merger. Various cuts (hacking attempts, browser 
misconfigurations, etc.) removed about 5% of our data. 
The data were reduced into two final catalogues based 
on whether data was weighted or unweighted. In the un- 
weighted data, each user's classification carries an equal 
weight, while in the weighted case, users weights are iter- 
atively adjusted according to how well each users agrees 
with the classifications of other users. In both cases, the 
accrued classifications are further distilled into super- 
clean, clean and cleanish catalogs of objects, for which we 
require 95%, 80% and 60% of users to agree on a given 
classification. In all cases, this is a statistically significant 
detection with respect to random voting; however, the hu- 
man "systematical" error associated with it is difficult to 
judge. In any case, we are in the limit where taking more 
data will not change our sample be yond noise fl uctua- 
tions as the votes are uncorrelated. In lLand et al. 

I i2008) 

a bias of unknown origin toward anti-clockwise galaxies 
was discovered and corrected for by adjusting the clean- 
liness level for the clock-wise galaxies to a slightly lower 
value. This work uses the same data and bias correc- 
tion. We note, however, that if unaccounted for, such bias 
would generate a constant offset in the correlation func- 
tion that cannot mimic the correlations we are seeing in 
the data. After bias correction is applied, the numbers of 
clockwise and anti-clockwise galaxies are the same within 
Poisson noise in each sample. 

We decided to use the 80% clean weighted sub- 
sample. We stress that the decision to work with 80% 
clean sample was made in advance and was not chosen 
to maximize our signal. We show an example of a typi- 
cal clockwise and anti-clockwise spinning galaxies from a 
clean catalogue in the Figure [1] 



4.2 3D, angular and projected configuration 
spaces 

In the formalism of Section we have always referred 
to distance between two galaxies as being r, the physical 
distance in real-space. In practice, it can be any measure 
of distance between galaxies. In this work we use three 
different distance measures: 

(i) Angular distances. These have the advantage of 
producing the highest number of pairs. We denote the 
corresponding correlation function with c{8). 

(ii) Real space distance. We use the distance in the 
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Figure 1 . This figure sliow a pair of typical galaxies from our 
clean catalogue. The left image is an anti-clockwise (S-like), 
while the right is a typical clockwise (Z-like) galaxy. 



redshift space for pairs of galaxies for which both spec- 
troscopic redshifts are known. These are not the true 3- 
dimensional distances, but are instead distances in the 
red-shift space an d therefore affec ted by the fingers-of- 
god effects (see e.g. lHamiltonl ll998V Since the axis of sub- 
halos are correlated with the sh ape of the parent halo (see 
e. g. iBailin fc Steinmetd[2005h . there exist correlations 
in the ratio of edge-on to face-on spirals as a function 
of projected distances from the centre of the halo. This 
considerably complicates any correction for fingers-of-god 
effects and therefore we do not attempt this correction, 
since effects are likely to be sudominant. A concordant 
flat cosmology with f2m = 0.25 was assumed when calcu- 
lating distances. We denote the corresponding correlation 
function with c(r). 

(iii) Projected distances. There distances are the 
transversal component of the distance vector connecting 
two galaxies with known redshift. If only one galaxy in 
the pair has a known redshift, we assume the other galaxy 
to have the same redshift. The advantage of this distance 
measure is that it is not affected by the redshift space 
distortions and that the number of pairs is signiflcantly 
larger than in the case of real space distances. We denote 
the corresponding correlation function with c{p). 

For each of the above distance measures, we flrst 
located all galaxy pairs in our sample, that are less than 
2000 arc seconds or 3 Mpc//i or 1 Mpc//i projected apart. 
This gave us three sets which we describe in the Table [T] 

We have then removed rogue pairs. In the primary 
SDSS pipeline analysis, every object is assigned an SDSS 
ID. Large nearby galaxies are often associated with more 
than one ID, as various knots and substructure of the 
galaxy are recognised as sources by the reduction soft- 
ware. All such IDs are therefore classified as the same 
galaxy, resulting in spurious positive correlation at the 
shortest distances. Our automatic mechanism removed 
all pairs for which their angular separation is less than 
1.5max(rp), w here max(rp) den otes the larger of the two 
Petrosian radii l|Petrosianll 19761 ) . This did remove the ma- 
jority, but not all of the rogue pairs. Therefore, the clos- 
est pairs (at angular separations of less than 3rp) in each 
category were examined by hand and 69 additional SDSS 
objects were removed. 



4.3 Determination of a angle 

As discussed in the Section O we need to estimate the 
value of a, the maximum angle of inclination at which 
the spirals have a measured spin orientation rather than 
a being classified as "edge-on" spirals. To do this, we use 
the adaptive second moments (|Bernstein fc Jarvisll200 
from the SDSS pipeline , namely ex a nd e+ to calculate 
the axis ratio, following iRvdenI (|2004| ): 



9 = 



1 - e 
1 + e 



1/2 



(25) 



where e 



+ . In Figure [2] we show the distribu- 
tion of q values for spirals galaxies classified as face-on (of 
either spin orientation) and edge-on spirals. As expected, 
the two populations occupy the two corners of possible 
values of q, but there is a significant overlap. Intrinsic 
ellipticities, non-zero thickness of the disk and potential 
human-induced selection effects likely complicate things. 
We have attempted to model in trinsic ellipticities in the 
spirit of iGiovanelli et al] l|l997h . but difference was neg- 
ligible. 

A plausible range of the cut-off q is 0.2 - 0.5, giving 
the values of a between 60° and 80°. If we numerically 
integrate Equation (|13[> as explained in Section 13.11 we 
get 



(26) 



with the value of m between ~ 0.6 and ~ 0.9. We will 
assume a systematic bias associated with this effect to 
be m = 0.75 ± 0.15. Adding to this the effect of the 
random reversing of galaxy spins and allowing a liberal 
50% enhancement of the systematic error due to an ad- 
hoc assumption in Equation (|15p . we arrive at 



r7(r) = (0.9±0.3)c'"'="=(r). (27) 



4.4 Correlation function measurement 

Our basic method is to measure c(r) and its errors and 
then to infer constraints on ri{r). 

To measure c(r), we note the following. For a pair 
of galaxies, whose spins are Si and Sj, the product SiSj 
can be either -1-1 or -1 with probabilities p±i- Since + 
P-i = 1 and the expectation value of (siSj) = p+i— p_i = 
c(r) it follows that 



P±i 



1 ± c(r) 



(28) 



Therefore, one can write the likelihood function for 

c(r) as 



P(c(r) |data) oc P(datalc(r)) = 



1 + dkc(rk) 



(29) 

where index k runs over all pairs of galaxies in the 
sample and dk = SiSj is the spin product for the fc-th 
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Property 



Angular 



Real space 



Projected 



Number of pairs 
Number of gal. 
Mean z 



34031 
20827 
0.08 



8005 
7979 
0.05 



24271 
25272 
0.07 



Ax^ exponential 
Evidence ratio 



0.94 



9.76 
9 

+0.39+0.54+0.56 
-0.48-0.80-0.92 



0.35 



9.12 

6 

+0.19+0.42+0.86 
0.16-0.27-0.34 



0.55 



5.89 
1 

+0.60+0.89+0.95 
-0.47-0.54-0.55 



23.41 



hll. 37+42. 47+73. 24 
-6.35-13.04-22.41 



0.37 



+0.16+0.41+0.61 
-0.11-0.23-0.37 



0.02 



+0.08+0.37+0.47 
0.01-0.02-0.02 



Ax^ Gaussian 
Evidence ratio 



0.60 



9.82 
11 

+0.25+0.38+0.40 
-0.27-0.49-0.58 



0.24 



11.52 
16 

+0.10+0.22+0.41 
-0.09-0.17-0.23 



0.44 



6.71 

1 

+0.31+0.51+0.55 
-0.37-0.43-0.44 



26.76 



+9.36+34.11+68.56 
6.35-13.71-25.42 



0.44 



+0.13+0.31+0.52 
0.10-0.20-0.42 



0.02 



+0.05+0.36+0.47 
0.01-0.02-0.02 



Table 1. This table shows the basic information about the datasets used in this work. When calculating the mean redshift, only 
subset of galaxies with redshift is used and we average over galaxies and not galaxy redshifts. We also report values of best fit 
Ax^, Bayesian evidences and parameters of our fits. Note that evidence here is the evidence ratio and not its logarithm. 



0.05 



0.04 - 




q 

Figure 2. This figure shows the histogram of distributions of 
q values for galaxies classified as face-on spirals (solid black, 
classifications 2,3) and edge-on spirals (dashed blue, classifi- 
cation 4) 

pair whose distance is r^. In practice we work with the 
log likelihoods 

logP(c(r-)|data) = ^log(l + djc{rj)) + const. (30) 

We use three possible forms for c(r). First we assume 
a stepwise shape for c(r) and measure it in bins. Second, 
we use two 2-parameter families of curves that seem to 
describe our data fairly well: an exponential 

c(r) = min 1 1, ae~'^^'j> (31) 

and a Gaussian 

c(r) ^ae-'-'/""'. (32) 



This parameter space is so small that it can be effi- 
ciently explored using grid based methods and more ad- 
vanced Markov Chain methods are not necessary. 



5 RESULTS 

In Figure [3] we plot the results of our binned estimation 
of c(r). From the two figures it is immediately clear that 
there is a hint of an excess at low values of r. The statis- 
tical significance of this excess is marginal, at about Ax^ 
of 7.5, 14.2 and 5.6 for angular, real and projected dis- 
tances with 6 extra degrees of freedom associated with 6 
bins. This corresponds to 2-3 cr detection in the redshift 
space but a non-detection in other spaces. 

To understand this excess better, we calculate the 
probability contours on the a-b plane using exponential 
and Gaussian likelihoods. These are plotted in the Figure 
|4]and the relevant numbers are in the Table [T] How sig- 
nificant are these detections? The improvement in is 
between 9 and 12 with respect to zero correlation in angu- 
lar and redshift cases with two free parameters. Within a 
frequentist approach this is significant at 2-3 sigma level. 
The excess at low redshift is not significant in the case 
of projected distances, although visually the low distance 
points are not incompatible with an excess. 

A more appr opriate statistical procedure is the 
Bayesian evid ence (j Slosar et al. 2003; Boltran ct al. 2005; 
iTrottal 120071 ') which we calculate for all our 2 model 
parametrizations and also show in table [T] These can 
be calculated exactly for a simple problem like ours. 
Evidence depends weakly on the prior size and in this 
we chose the prior on a between and 1/1.5 for Gaus- 
sian/exponential case and b between and 1000 arc sec 
or 1 Mpc//i or 0.5 Mpc//i projected. Regardless of the 
exact number employed, the evidence ratio is between a 
few and a few tens units implying a weak evidence or 
a hint for angular and redshift spaces, but not for the 
projected space. This is consistent with results from the 
frequentist approach above. 
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Figure 3. This figure shows the constrains on the binned 
correlation function c for angular (top), redshift (middle) and 
projected (bottom) spaces. Two lines correspond to our best 
fit exponential (solid red) and Gaussian (dashed green) fits. 



Finally, we acknowledge the fact that the exponential 
and Gaussian form were chosen a-posteriori, after seeing 
the data and hence the improvements in fits contain a 
subjective a-posteriori factor. 



5.1 Systematics 

We can now briefly discuss some of the main systematic 
effect that might affect our measurements. 

Rogue pairs. As discussed in Section 14.11 we manu- 
ally looked at all pairs in the clean sample and discarded 
rogue pairs. It is an important systematic check, because 
we have at the same time convinced ourselves that man- 
ually classifying a small subset (80 galaxies) of the total 
sample gave consistent results. 

Weighting. Repeating our measurements with un- 
weighted data, changes results by less than 5%. 

Cleanliness level. We have repeated the analysis with 
the super-clean sample. There are many fewe r galaxies in 
the super-clean sample (|Lintott et al.ll2008l ) and so the 
statistical signiflcance decreases considerably. We have no 
significant detection in any of the spaces considered. The 
errorbars increase by a factor of 2 to 2.5, but the central 
values in individual bins remain consistent. While the 
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Figure 4. This figure shows the constrains on the a-b plane 
for all datasets and models under consideration. Thick lines 
enclose 68.3%, 99.4% and 99.7% enclose likelihood contours for 
the weighted sample. Thin lines are the same for unweighted 
sample. The top and bottom rows show results in real and 
angular spaces respectively. The left and right columns the 
exponential and Gaussian fittings exponentially. 

statistical power is decreased, the final signal is consistent 
with the results presented above. 

We have also repeated our measurements with the 
the cleanish sample that requires 60% of votes to agree. 
The results imply strong detections in both angular and 
projected samples, but with a lower significance in red- 
shift space. We show their results in the Figure [S] The 
high-confidence with which the results are detected in 
angular and projected spaces is likely to be deceiving, as 
the rogue pairs have not been manually cleaned for these 
samples. The decrease of signal in the redshift space indi- 
cates that the signal is indeed getting lower due to noise 
introduced by low significance. 

Selection effects. Another important question is 
whether there is any physical difference between our red- 
shift sample and angular sample and how do these sam- 
ples compare to the general SDSS sample. To do this 
we have divided the galaxies that formed our pairs at 
closest distances into those for which we have redshift 
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Figure 5. Same as Figure |3] but for cleanish sample. These 
results are likely to be affected by the rogue pairs - see text 
for discussion. 



information and those for which redshift information is 
not available. When comparing colors and magnitudes 
we find that there is no evidence that objects with and 
without redshifts are drawn from different magnitude, 
u — r colour or petro radius size distributions. The rea- 
son for some objects lacking redshift is therefore probably 
incompleteness due to fibre collisions. We therefore find 
no evidence that the correlation in the angular sample is 
of a different physical origin than the correlation in the 
real-space sample. 



6 DISCUSSION AND CONCLUSIONS 

We are now in position to make a synthesis of our results. 
The redshift-space results show that there is a significant 
correlation of c(r) ~ 0.15 in the projected galaxy spins 
up to the ~ 0.5 /i/Mpc. The angular correlation, shows 
larger correlations of c{9) ~ 0.4 that are significant corre- 
lations up to 30 arc-seconds, which roughly corresponds 
to projected distances of about 0.03 Mpc/h, since mean 
redshift of 0.08 correspond to distance of ~ 230 Mpc/h. 
This is consistent with the redshift-sample results - both 
exponential and Gaussian fits do predict c(r) to raise to ~ 
0.3-0.4 as r goes to zero. A consistent picture is therefore 
the following: The angular sample detect correlations at 



the shortest distances, where majority of pairs are physi- 
cal associations, but these get diluted at larger distances 
due to interlopers. The redshift-space correlations track 
these correlations to larger physical distances. Projected 
space pairs do not have enough signal-to-noise to detect 
these correlations at high significance. 

How do these results compare to theoreti cal predic- 
tions? Simple models as those suggested in I Pen et all 
|2000l ) (equations Q and ([5])) predict a vanishing r; and 
hence we have directly detected a non-random distribu- 
tion of inertia tensors. Within the standard model, the 
reason for correlations of moments of intertia are the cor- 
relations of these with the (slowly varying) tidal field. 
On the other hand, if moments of intertia are perfectly 
aligned with the tidal field, the tidal torque cannot pro- 
duce any angular momentum and therefore the resulting 
angula r momentum is due to the residual 10% of misalig- 
ment (|Porciani et al.l l2002lj) . The stunning outcome of 
our result, if confirmed, is that even these 10% misalig- 
ments are correlated from (sub-)halo to (sub-)halo. 

What is also i nteresting is, however, that in 
IPorciani et al.l l|2002al ) 77 correlations have not been de- 
tected in simulations at 2: = at all separations. In par- 
ticular, ri < 0.02 at r = 0.5 Mpc/fe. A virtually ide ntical 
result has been found bv lBailin fc Steinmetd ()2005h . who 
also find rj < 0.02 at r = 0.5Mpc//i (our 77 is their 
This is in tension with our results even after conversion 
factors in Equation (|27|l are taken into account. There 
are many reasons that explain why our results are not 
directly comparable to the above work. First, they are 
comparing individual dark matter halos. In our case, we 
see the signal at pair separations of less than 1 Mpc. At 
such distances, one-halo pairs (pairs of galaxies that re- 
side in the same dark matter halo), dominate over two 
halo pairs (pairs in which two galaxies occupy two dif- 
ferent halos). By selecting spiral galaxies, we are essen- 
tially selecting pairs that are composed of satellites re- 
siding in the same halo, rather than pairs compromised 
of central halo galaxies. The latter are bright ellipticals 
and hence inaccessible using our method. Unfortunately, 
not very much theoretical work has been done for spin 
correlations of s ubstructur e. The most relevant paper in 
the literature is iLee et all (|2005l ). which, however, still 
uses the chirality agnostic model of lPen et all l|2000l ) and 
does not calculate the chiral correlation function. More 
work on the theoretical side and A'-body side is required 
to understand the implications of our results. Hopefully, 
the results could be turned around and help us under- 
stand what kind of substructure spiral galaxies occupy in 
a typical dark matter halo. 

It is therefore imperative that our observational pro- 
tocol is simulated on a large enough A-bo dy simulation, 
for e xample the Millennium Simula tion JSpringel et al 



'2005') or MareNostrum Universe (|Gottlober fc Yepei 



;2007. ) simulations. There, halos and sub-halos hosting spi- 
ral galaxies can be identified and those, whose inclination 
with respect to a given observer is small enough to be 
considered face-on, should be correlated. This would re- 
sult in a quantity c(r) that is directly comparable to the 
observables that we constrain with the Galaxy Zoo data. 

Another interesting aspect of our results is that, for 
spiral galaxies, we essentially exclude large and random 
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misaligments between gas and dark matter angular mo- 
menta. Since the dark matter is dynamically dominant, 
gas angular momenta can only be correlated if they are 
so due to correlations between dark matter. 

Finally, it is tempting to combine our measurements 
with the ellipticity measurements to improve signal-to- 
noise and remove some systematic. Note, however, that 
our 1-bit signal divides a four-fold degeneracy into a two- 
fold one and thus this is a non-trivial task which will be 
left for the future. 

To conclude: We have tentatively detected a chi- 
ral correlation function in the spins of spiral galaxies. 
This correlation function vanishes in the simplest models 
based on tidal torque theory. Our results indicate, that 
moments of intertia of protohalos that end up hosting 
spiral galaxies are correlated at distances less than ~ 0.5 
Mpc//i. These short distances imply that these protoha- 
los are often likely to be substructures of massive halos. 
More work is required to understand these results at a 
quantitative level. 



ACKNOWLEDGEMENTS 

We would be nowhere without the amazing contributions 
from all the GZ members, forum contributors, and other 
essential volunteers! 

AS thanks Cristiano Porciani for bringing up the im- 
portant distinction between one and two-halo pairs and 
Uros Seljak for useful literature tips. 

AS is supported by the inaugural BCCP Fellowship. 
CJL acknowledges support from an STFC Science in So- 
ciety Fellowship. 



APPENDIX A: CALCULATIONS OF "SEEN 
PAIR" PROBABILITIES. 

The question that we want answer is; For a given pair 
of galaxies, whose spin vectors are at angle cos 6 = /i, 
what is the probability of an observer seeing the pair 
with the same sense of galaxy rotation, what opposite 
senses of galaxy rotation and not seeing them at all due 
to selection effects? 

If Q = tt/2, the result can be obtained by considering 
each spin in turn. The first spin divides the unit sphere 
of possible observer directions into two half-spheres, de- 
pending on the sign of its projected spin. When two spins 
are considered, the intersection of the two half-spheres are 
two lunes. The thickness of the lune of opposite spins is 
9/tt, leading to the result in Equations and (|12p . 

When a < n/2 the dividing line between the two 
half-sphere becomes a band of angular thickness 2{n/2 — 
a) and the two half-spheres shrink to two spherical caps 
of radius a. The overlapping area of the spherical caps 
separ ated by 6 is 47rf(0), w here f{9) is given by Equation 
([2^ (|Oat fc Sandeill2007l ). 

If a > 7r/4, the intersection of the 4 cups give four 
"trimmed" lunes. There are three possibilities. 

• S < 2(7r/2 — q). Both spins are in roughly same direc- 
tion and the opposite spin "trimmed" lunes are squeezed 
to zero area. Hence P+i = f{9) and P_i = 



• 2(7r/2 — a) < 6 < 2a. General situation in which all 
four "trimmed" lunes have finite area. We have P+i = 
fie) and P_i = fiTT - 6) 

• > 2a. Both spins are in roughly opposite directions 
and opposite spin "trimmed" lunes are squeezed to zero 
area. In this caseP+i = and P_i — /(tt — 6^) 

These results imply Equations pO)) and (f2T|) . We 
have tested these analytical predictions using a Monte 
Carlo code. 
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